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SUMMARY 


A  new  set  of  phase  variables  is  proposed  and  justified  to 
develop  automatically  by  computer  the  solar  terms  in  Lunar  Theory. 

The  dependence  of  the  perturbation  function  on  the  mass  ratio 
Moon/ (Earth  +  Moon)  is  completely  elucidated.  A  recursive  procedure 
is  proposed  to  develop  that  function  so  as  to  keep  explicit  all 
its  d'Alembert  characteristics. 

The  perturbation  series  obtained  by  computer  is  compared  with 


Delaunay's  development. 


By  an  Analytical  Lunar  Ephemeris,  we  mean  a  general  approximate 
solution  applying  to  any  single  moon  whose  orbit  satisfies  certain 
a  priori  assumptions;  the  coordinates  and  velocity  components,  or  any 
other  state  variable,  come  out  as  literal  expressions  in  terms  of 
initial  constants  and  of  dynamical  parameters  carried  throughout  the 
theory  as  symbols  to  which  nowhere  is  given  a  numerical  value.  An 
Analytical  Lunar  Ephemeris  presents  definite  advantages.  The  degree 
of  accuracy  to  which  the  theory  is  constructed  is  given  an  explicit 
definition  in  reference  to  the  physical  assumptions  and  the  size  of 
the  parameters  involved.  It  is  immediately  applicable  to  any  single 
moon  in  our  solar  system.  It  is  arranged  in  such  a  way  that  any  small 
change  caused  hy  improved  data  on  the  values  of  the  constants  can  be 
made  easily  without  requiring  to  go  over  the  whole  of  the  work  again. 

In  the  operations  of  fitting  the  theory  to  the  observations ,  it  readily 
provides  the  true  partial  differential  coefficients  relative  to  the 
model  on  which  it  was  constructed.  It  is  the  only  sure  way  to  decide 
whether  the  residuals  0-C  lie  in  inadequacies  of  the  nominal  model  or 
in  inaccuracies  of  the  solution.  Finally  a  numerical  or  operational 
mistake  in  the  course  of  building  an  algebraic  expression  is  very  easily 
traceable. 

On  the  other  hand,  the  number  of  terms  which  have  been  found 
necessary  to  secure  a  degree  of  accuracy  commensurate  with  that  of 
observations  is  very  large,  and,  without  the  help  of  a  computer,  it 
becomes  an  impossible  task  to  obtain  them  with  any  degree  of  certainty 

and  thoroughness.  PRECEDING  m  BLANK 


The  formation  of  an  Analytical  Lunar  Lphemeris  may  be  roughly 
divided  into  three  stages.  At  first,  the  three  bodies  Sun,  Earth  and 
Moon  are  taken  as  mass  particles;  the  barycenter  of  the  pair  Earth- 
Moon  is  supposed  to  move  round  the  Sun  in  a  given  Keplerian  orbit;  all 
displacements  of  this  orbit  and  of  the  Moon's  orbit  caused  by  any 
other  source  than  the  Sun  and  the  Earth  as  mass  particles  are  ignored. 

This  Initial  stage  is  called  the  Main  Problem,  and  its  solution 
constitutes  the  Solar  Part  of  an  Analytical  Lunar  Lphemeris.  In  the 
second  phase,  secondary  perturbations  are  investigated.  The  algebraic 
construction  may  proceed,  at  least  in  principle,  without  any  knowledge 
of  initial  conditions  beyond  a  crude  estimate  of  the  relative  magnitude 
among  some  of  the  constants  involved.  The  third  stage  would  then  be 
devoted  to  mounting  analytically  the  partial  derivatives  witli  respect 
to  the  astronomical  constants  to  be  determined  from  observations,  in 
evaluating  the  constants  numerically  and  in  substituting  these  numbers 
in  the  analytical  expressions  so  as  to  obtain  a  Numerical  Lunar 
Ephemeris. 

Of  the  methods  which  have  been  proposed  to  solve  analytically  the  Main 
Problem,  Delaunay's  T)i6otne  du  Mouvement  da  la  Lune  must  take  the  first 
place.  His  algorithms  may  not  be  the  best  adapted  to  automated  algebraic 
calculation  (Barton  1966),  but  he  has  actually  carried  his  construction 
to  a  high  degree  and  with  a  detail  greater  than  any  other.  Delaunay's 
formulae  offer  numerous  opportunities  for  testing  subroutines  and  assess¬ 
ing  what  new  analytical  solutions  of  the  Main  Problem  may  achieve. 
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During  the  last  four  years,  we  have  been  considering  implementing 
by  computer  the  suggestion  of  Dr.  Brouwer  (1961).  The  algorithm  of 
Von  Zeipel  which  h*»  proposed  may  appear  troublesome,  chiefly  from  its 
inability  to  produce  the  eliminating  canonical  transformations  and 
their  inverses  in  an  explicit  form.  Mersman  (1969)  has  significantly 
developed  the  capabilities  of  Von  Zeipel's  method.  Yet  Mersman's 
algorithms  for  unmixing  the  variables  are  not  recursive,  which  makes  it 
difficult  to  involve  them  in  general  codes  where  the  order  of  a  theory 
enters  as  dummy  argument.  Some  time  ago  (Deprit  1969),  it  has  been 
shown  how  the  difficulties  raised  by  Von  Zeipel's  method  might  be 
avoided  by  the  use  of  hie  transforms.  At  the  same  time,  our  computer 
techniques  (I)anby  et  al.  1965,  Rom  1970)  have  undergone  a  complete 
overhauling  to  accommodate  divisors  made  of  products  of  linear 
combinations  of  mean  motions. 

It  is  intended  to  combine  these  two  techniques — the  canonical 
perturbation  theory  based  on  Lie  transforms  and  the  processing  of 
bcheloned  Series — so  as  to  completely  solve  the  Main  Problem  of  an 
Analytical  Lunar  bmphemeris. 
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* 


1.  THE  PERTURBATION  FUNCTION 

Let  (wep.  j^,^)  designate  the  position  of  the  Earth 

(wop.  of  the  Moon,  of  the  Sun);  to  these  position  vectors  are 
associated  the  impulsions 

Jii  ■  (i-0,1,2)  (1) 

in  which  mQ  (wap.  m^n^)  indicates  th»  mass  of  the  Earth 
(wap.  of  the  Moon,  of  the  Sun).  The  mutual  distances  being  the 
functions 


ro,i "  ll^ill*  ro,2 "  ll-vM’  ri,2  *  (2) 


this  problem  of  three  bodies  is  described  by  the  Hamiltonian 

1  ii..  1 1  2  .2 


l  m 

2  0<i<2  1 


jiiir-k4  ^  -j-1 

(i,j)  ij 


(3) 


Positions  (u^,^,^)  in  the  inertial  frame  are  transferred  to 
positions  in  a  chain  of  barycentric  coordinate  systems 

according  to  the  formulas 


<&1  "  ~  ~0’ 

2S2  "  H2  ~  (1_al)'ii0  ~  °1*1»  (4) 

2Sq  *  a2«y<2  +  Ci-a2)  [  (i-ai)^  + 

whe rein 

ni^  ^2 

^1  mQ+m^  ’  °2  mQ+m^+n^  '  ^ 


The  inverse  mapping 
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*0  m&0  "  °1*1  "  0  2&2  * 

^1  -  &0  +  “  °&2*  <6) 
U2m&0+  (1“02)<^2 

extends  through  the  implicit  definitions 


J!o  -  a-opa-oj)^  -  *{  -  al4j, 


%  -  <’1a-»2tt0  +asf  -  u-oj)^. 


(7) 


jj2  -  “A+JS? 

into  a  homogeneous  linear  canonical  transformation.  The  vector  ^ 
is  the  geocentric  position  of  the  Moon,  the  vector  jt.,  the  position 
of  the  Sun  with  respect  to  the  barycenter  of  the  couple  (Earth-Moon) 
and  the  vector  jCg,  the  inertial  position  of  the  barycenter  of  the 
system  (Earth,  Moon,  Sun),  Under  this  homogeneous  canonical  extension 
the  Hamiltonian  decomposes  into  the  sum 


mnmi  ("U+m.  )m_ 

»  -  en  +  +  3»* 

0  nin+in.  1  mri+m,+m„  2 


0  1 


0  1  2 


(8) 


whose  terms  are 


CQ  -  |(m0+m1+m2)||  Xq  ||2, 


"Til 


m  +m  k  (mn+m. ) 

x*  I!2 - 2_i_ 

m0ml  1  "  rl 


<V 


(92) 


m0+m1+m2  2  k  (mg+m^+m,,) 

(m0+m1)m2  ^2  II  r2 


(93) 


T*  =  k2m0[m~(— - - — )  +  m.  (— - - — 

2  °  r2  r0,2  1  r2  r0,l 


)]. 


(V 


The  distance  functions  entering  the  Hamiltonian  are  now 
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From  now  on  the  orbital  plane  of  the  Sun  is  taken  as  the  reference 
plane  (x,y).  If  g2  is  the  argument  of  perigee  of  the  Sun  reckoned 
from  the  x-axis  and  f2  its  true  anomaly,  then  the  coordinates  of 
the  Sun's  position  £2  are 


r2cos(f2  +  g2). 

<15^ 

r2sin(f2  +  g2)» 

(152) 

0. 

(153) 

Under  this  approximation,  the  motion  of  the  Moon  is  described  by 
the  Hamiltonian 


M  -  *,£1  J&2  ^  +  ^i**2(t)) 


with 


2  k  (m()+m1) 


»-  (7-  -  71-)  +  77-  (7-  -  t*-)]. 

num f  4  o j  1*2  Lq  2  * *0 2  r 2  » |  2 


“0  1 


(16^ 


(162) 


(163) 


But,  in  application  of  the  binomial  law  under  the  assumption  that 


rl  <  r2' 


— —  I  <-o,)n(— )*?  (s), 

r0,2  r2  „>0  1  r2  " 


-i-.ij  <  >«£)%  <s> 

rl,2  r2  n>0  1  r2  n 


where,  for  any  n  >_  0,  is  Legendre's  polynomial  of  degree  n. 


Consequently 


k2m„ 


3>=-  — -  v  [(l-o  )n_1  -  (-a  )n_1]A)np  (s).  (17) 

r2  nV2  1  1  r2  n 


2.  BASIC  PHASE  COORDINATES 

We  substitute  Delaunay’s  phase  variables 
to  the  Cartesian  coordinates  £  ■  (x^.y^.z^)  and  conjugate  momenta 
■  (Xj,Yj,Zj).  Then  we  modify  them  into  the  set 


pi  '  ‘Ei  •  V 

8i "  "hr 


(18) 


The  phase  space  of  the  Main  Problem  in  Lunar  Theory  will  be  here 
coordinatlzed  by  this  set  of  modified  Delaunay's  elements. 


In  relation  with  Delaunay's  action 
axis  a^  and  the  mean  motion  n^  such  that 


2  2 

Aj_  ■  k  (mQ  +  m1)a1, 
2  3  2 

i^a,  -  k  (mQ  +  n^). 


we  define  the  semi-major 

<19^ 

(192) 


Then  we  Introduce  the  parameter  of  lunar  parallax 


the  modified  parameter  of  lunar  parallax 

*i  ■  (1_2° i> j ; 

also  we  designate  by 

m  *  n2/n^ 


(20) 


(21) 


(22) 


the  ratio  of  the  mean  motions  of  Sun  and  Moon. 
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In  relation  with  the  norm  of  the  angular  momentum  G^,  we  found 
convenient  to  use,  beside  the  eccentricity  e^,  a  function 
E^  ■:  E^(L^,P^)  >  0  such  that 

2PX  -  I^E2.  (23) 

From  the  obvious  relation 


E2  .  I  E« 

1  4  Ll* 


there  results  that 


IE2 

2  L1 


1  -  <l-e2)1/2. 


(24) 


(25) 


In  terms  of  what  is  sometimes  called  the  anomaly  of  the  eccentricity, 
namely  the  angle  such  that 

0  <  ^  <_n/2,  *  sin  4^, 

it  turns  out  that 


E^  ■  2  sin(4>^/2) . 

In  relation  with  the  component  of  the  angular  momentum 

normal  to  the  ecliptic,  we  found  convenient  to  use,  beside  the 
inclination  1^,  a  function  _  J^(L^,Q^)  >  0  such  that 

2Q1  -  I^J2. 

It  is  easily  seen  that 

1  o/i  2*1/4  .  1  t  _  1  r2*l/2  ,  T 

m  2(1  e^)  sin  2  I 2(1  ~  2"  E^)  sin  ^  I^» 


(26) 
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Also,  we  shall  often  substitute  to  the  angles  g^  and  h^  the 
linear  combinations 


F1  “  £1  +  81»  (271} 
D1  "  X1  ‘  X2  “  ll  +  81  +  hl  “  *2  "  g2  (272} 


that  have  been  introduced  by  Delaunay  under  the  names  of  mean  elonga¬ 
tion  of  the  Moon  to  its  node  and  of  mean  elongation  of  the  Moon  to  the 
Sun,  respectively. 


The  development  of  the  Analytical  Lunar  Ephemeris  involves  a 
number  of  basic  differentiations.  Reviewing  them  here  will  justify 
partly  why  we  departed  from  Delaunay's  traditional  functions  e^  and 
1^  in  favor  of  the  quantities  and  J  . 

From  the  elementary  derivatives 


_L_ 

3L. 


nlSl 


3 

3L, 


» 


we  immediately  deduce  that 


k  (mg+ti^) 


and,  for  any  j  >  0, 


3L1  al 

iv 

3  1  3 

3L^  n^ 

- 

o, 

2 

k  (m0+m1)  , 

3 

a.  3L. 

1 

1 

2 

k  (m0+m1)  , 

3_ 

-2n^ , 


or 


-l 


a*  *  2jn1oi^, 


J. 


9LX  m  =  3jnlal; 


(281) 

(282) 


(283) 

(284) 
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(29x) 

(292) 

<301) 

(302) 


We  conclude  this  sequence  of  elementary  differentiations  by  a  set  of 
simple  formulas  for  the  partial  derivatives  with  respect  to  the  actions 
Ap  ^1’  ^1  t*ie  two  tyries  monomials  that  we  shall  meet  time  and 
again  in  the  development  of  an  Analytical  Lunar  Ephemeris. 


(nu+m 
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3.  DEVELOPMENT  OF  THE  PERTURBATION 

Let  denote  the  true  anomaly  of  the  Moon;  then  its  position 

vector  has  for  components 

■  r1 [cos (f ^fg^Jcos  h,  -  sin^+g^sin  h^  cos  1^],  ( 32 ^ ) 

y x  a  r1[cos(f1+g1)sin  +  sin (f^+g^cos  h^  cos  1^],  (32 

z1  =  r1sin(f1+g1)sin  1^  (323> 

and  the  elongation  function  s  defined  in  (11)  turns  out  to  be 
the  function 


s  =  cos[(f1~£1)  -  (f2'^  +  Di^cos2  2  *1 

+  cosKfj-f.j)  +  (f2-i2)  +  2Fj_  -  D1]sin2  ~  I  . 

(33) 

The  systematic  development  of  s  may  proceed  as  f ollov;s  . 

(i)  One  obtains  the  d'Alembert  series 


cos (f-;) 


1  +  N 


cos 


sin(f-£)  = 


1  (  1  sa  ke2k\ejsin  j£. 
jil\k>0  j 


(34x) 

(342) 


Several  algorithms  have  already  been  proposed  to  generate  such  series 
by  computer  (Deprit  and  Rom  1967,  Broucke  1969);  we  propose  to  come 
back  to  this  elementary  but  important  problem  in  the  next  installment. 
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(11)  In  the  series  (34)  applied  to  the  Moon,  i.e.  for  f  -  f^, 


£  ■  £^,  and  e  -  e^,  we  make  the  substitutions 


2n  /t,2  1  r,4Nn 


el  "  <E1  “  f  El>  * 


e2n+1  ■  (e2  -  —  E4)en 
el  U1  4  V  1 


where  is  the  expansion  of 


ei  -  Ei(1 


I  e2,1/2 

b  V 


in  powers  of  E^.  These  substitutions  produce  a  d’Alembert  series  of 


the  type 


<frV  * 1  +  cj,kEik)EiC0S  JV 


sin(f1-£1) 


2  ( I  si ,kEi  ) Eisin  J'r 

jnykio  / 


(iii)  Then,  by  straightforward  multiplication,  we  calculate  the 


products 


cos (f 1~£1)cos (f 2— ^2^ » 
cos (f^-£^)sin(f2-£2) , 


-  sin(f1~£1)cos (f 2~£2 ) , 


sin(f1-£1)sin(f2~£2) , 


which  come  out  as  d'Alembert  series  in  the  pairs  (E.^,2^)  and  (e2,£2). 
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\ 

(iv)  As  the  last  intermediate  step,  we  produce 


cos[  (fj-J!^)  -  +  Dl^ 

cos[(f1-H1)  +  (f2-£2)  +  2?1  - 


(Mj+M^cos  Dx  +  (M2-M3)sin 
(M1-M4)cos(2F1-D1)  -  (M2+M3)sin(2F1-D1) . 


The  first  trigonometric  function  will  come  out  of  the  straightforward 
multiplication  as  a  series 


I 


h-0 


I 


i2lO 


X  I 
h  J2 


J 1 » J 2  ^2 

ci  i  E1  e2  C0S(J1£1  +  j2£2  -  V 


with  the  d'Alembert  characteristics 


jj_  5  (mod.  2),  | j1|  <_  i3, 

J2  -  i2  (mod*  2)*  |j2l  1  i2  5 


similarly  the  second  product  will  be  a  series 


12  2  1  c  1 

hi0  Si0  Ji  h  r 


J 1 » J 2  *1  *2 

i  E1  e2  COs(jl£l  +  j2l2  -  2F1  +  °1) 


with  the  same  d'Alembert  characteristics. 


According  to  (34),  the  elongation  function  s  can  now  be  assembled 
as  a  series  in  E^,  e2;  by  means  of  the  development 


,  2  1  t  1  .2 1  _2.-l  1  t2 

sin  —  I,  =  7  J,  1  -  r  E.  =v 

2  14  1  2  1  4 


A  2 

1  j>0  2J 


(35) 


it  can  be  reworked  as  the  series 
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^’1,^2,*3  ,^1  ^22^3 


II  I  II  I  c  T  i  E i  e 2  J i  cos(j  l  +j  i  +2j  F  +D  ) 

i^O  i2>0  i3-0,l  jx  j2  j3-0,l  J1’J2°3  111 


with  the  d'Alembert  characteristics 


I Jxl  1  ij_.  Jx  H  (mod.  2), 

I  Jo  I  1  V  Jo  -  io  (™>d.  2)» 


I J 3 1  —  ^3* 

At  this  point  we  depart  from  the  suggestions  made  by  Delaunay 
(1860,  p.  17  to  20):  a  systematic  procedure  is  in  order  if  the  development 
of  the  perturbation  is  to  be  carried  by  computer. 

From  Kepler's  third  law  applied  to  Sun  and  Moon,  we  have  that 


.2  2 
k  ra2al 


k  (niQ+n^)  2 

- m  a  . 

ai  2 


Delaunay  (1860,  p.  20)  decides  at  this  point  to  take  o 2  =  1.  Actually 
since  the  algebraic  manipulations  are  carried  by  computer,  there  is  no  need 
for  erasing  the  symbol  o2;  its  presence  in  the  calculations  is  very  useful 
for  preliminary  checks  and  for  tracking  the  origin  of  the  terms  appearing 
in  an  operation.  We  also  put 


rl  =  al°l’ 


r2  =  a2p2  ’ 

(l-o1)n“1-(-o1)“’1 

*n  (l-2a1)n_2 


for  n  >  2. 


(381) 


(382) 


(383) 
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With  these  conventions,  the  perturbation  function  3>  given  by  (17) 
can  be  slightly  reworked  to  appear  in  the  form 

/m  I  m  \  n  ^ 


k  (m^) 

* - -  °2m 


I  Vi 

J12  2 


1-2  1 

i  !FtV8)- 


In  this  manner,  3>  is  more  akin  to  de  Pontecoulant *9  perturbation 
than  to  Delaunay's  series:  the  parallactic  parameter  does  not 

enter  explicitly  the  part  of  *P  corresponding  to  j  =*  2  (Brown  1896, 
pp.  7  and  82). 

Let  us  examine  the  mass  functions  ;  to  this  effect,  we  introduce 


the  correction  factor 


o.U-Oj) 

,  -  — . .  ■  '  SS . -  ■  ,  ■  t 

1  (1-20J)2  (m0-m1) 2 


By  putting  also 


*  (1+4^)' 


and  inverting  the  definition  (40) ,  we  recognize  easily  that 


al  2B. 


1-0 l  '  26, 


X-2°l  *  t 


which  implies  that,  for  j  _  2, 


' ,  — 3ZT  1  +  <-l>X.  I 

1  2-1  0^k<J-2  K 1  1 

Since  they  constitute  polynomials  in  the  even  powers  of  the 

mass  functions  *  are  polynomials  in  the  correction  factor  c^. 
The  first  eight  ones  are  as  follows: 


* 


A 
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%  ■  1  +  c!» 

<J>3  *  1  +  2c^, 

2 

<t>6  =  1  +  3c1  +  c^, 

*  =  1  +  4ci  +  3c*, 

4>8  =  1  +  5c!  +  6cJ  +  cj, 

4>  *»  1  +  6c^  +  10c^  +  4c^. 


The  tradition  in  Analytical  Lunar  Theories  has  been  to  drop  from  CP 
the  terms  having  as  a  factor  (Brown  1896,  p.  8).  For  the  time 

being  we  shall  refrain  from  doing  so. 

Let  us  now  establish  a  recursive  procedure  to  expand  the  terms 

^ -4?rV*>-  <3-0)  <4 

From  a  library  of  elliptic  expansions  In  the  problem  of  two  bodies, 
we  produce  the  three  d  Alembert  series 


CP.  -  7-  -  1  +  I  I 


0  p 


|,(4 


»  k  •  I  ■;>  *  !,(!„  “>• 


il,i2r.il  i2. 


2  2  2  2  c'Ve^-dl1!^^- 

p2  ix>0  i2>0  jj_  J2 


(441) 

(442) 

(443) 


hence,  by  straightforward  multiplication  of  (36),  (442)  and  (443) 


the  d'Alembert  series 


r 


-19- 


11  1  II 

i£0  i^O  i3-Otl  j:  J2 


V  iri2,13pil  i2.2i 
j3  CJi*h,J3  1  2  1 


3 


costjj^  +  J2£2  +  2 J  3Fi  t  Di)‘ 


Remembering  the  recursive  relation 


(J+1)Pj+l(s)  "  (2j-,-1)8Pj(8)  “  JPj-i*  <J  1 
between  Legendre's  polynomials,  we  infer  the  recurrence 

(j+l)3>  (s)  -  ^  [(2j+l)s2>  -  j  ~  V).  (J  >_  1) 

j+i  P2  J  P2  J  1 


It  is  now  a  matter  of  straightforward  multiplications  and  linear 

combinations  of  Poisson  series  to  obtain  the  functions  in  the 

n 

form 


I 

i2l° 


i,  i~ 
I  E,2e,3 


i,>0 
^ _ 


1  2 


I 


I  1 

^2  ^3 


n>^2’^3*^4 

2  c 

Ja 


cos  (jjD^  +  +  J3i2  +  j^F^) 


with  the  d'Alembert  characteristics 


J1 5  n 

(mod .  2 ) , 

UJ 

1 n* 

J2  =  *2 

(mod.  2), 

1 J2I 

1  V 

^3  E  i3 

(mod.  2), 

U3I 

U  -  i4 

(mod.  2), 

1 J4I 

—  *4  ■ 

Substitution  of  these  terms  into  (39)  delivers  the  perturbation 

2 

k  (nig-hn^) 


T  -  T(a^t E^,J^,e2; 


,m,o2 ,c^ :D^,F^ ,1 j , S £ ) 
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as  .1  series  In  tin*  powers  of  the  modified  parallax  parameter  the 

Moon's  eccentricity  parameter  K^,  Its  Inclination  parameter  and 
tin*  Sun's  eccentricity  e 2 •  The  coefficients  in  this  series  are  finite 
trigonometric  sums  In  the  arguments  and  ^ *  the  pairs 

(•jfl’j).  <  !•  i  •  *  i )  *  ^u2,'2^*  present  the  d'Alembert  characteristics 

listed  for  any  J >  .  The  coefficients  of  the  trigonometric  sums  are 
polynomials  in  the  mass  correction  c^,  each  of  them  to  be  multiplied 
by  the  common  factor 

2 

k  (mQ+rn^)  2 
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4.  COMPARISON  WITH  DELAUNAY 

The  development  of  the  perturbation  has  been  performed  automatically 
on  an  IBM  360.  The  output  Is  a  deck  of  punched  cards  in  hexadecimal  code 
formated  to  serve  as  Input  in  the  further  construction  of  the  Analytical 
Ephemeris. 

The  main  program  and  its  subroutines  have  been  tested  by  reconstruct¬ 
ing  Jelaunay's  expansion.  Thus,  as  it  is  done  in  the  Tb'ovie  Ju  Vouvement 
>l<:  la  IjUuc,  we  eliminated  the  mass  ratios  by  assuming  that 

°2  "  *  °»  ci  "  0; 

also  we  used  the  eccentricity  e^  of  the  Moon  instead  of  the  function  E^, 
and  the  variable 

■i  ■  sIn  1  ‘i 

Instead  of  the  function  J^. 

Comparison  term  by  term  between  Delaunay's  text  and  the  computer 
listings  has  been  made  by  Mrs.  Deprit-Bartholome. 

In  the  expansion  we  had  retained  all  terms  whose  characteristic 

h  h  J3  h 
'  61  Y1  ®2 

is  such  that 

2 J i  +  +  J3  *  J4  -  10* 

Delaunay  has  explained  why  lie  did  not  collect  all  of  these  terms.  It 
may  be  of  interest  to  learn  by  how  much  Delaunay's  considerations  on  the 
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i 


"order"  of  vt)  have  simplified  his  perturbation  function.  We  enter 

2 

In  Table  I  the  103  missing  terms  having  <t  as  a  factor  (they  all  have 

characteristics  sucli  that  “  6)  and  in  Table  II  the  19 

3 

missing  terms  having  <  as  a  factor  (their  characteristics  are  such 
that  j9  +  =  4  and  J3  “  0)*  The  trigonometric  arguments  are  listed 

in  the  manner  Delaunay  ordered  them,  except  for  the  last  29  entries  in 
Table  I  and  the  last  8  ones  in  Table  II  which  are  not  to  be  found  in 
Del  tunay's  perturbation.  Some  of  the  terms  appear  in  Tables  I  and  II 
with  significantly  large  coefficients;  one  may  wonder  if  Delaunay  was 
justified  in  omitting  them. 

All  other  terms  that  are  common  to  Delaunay's  perturbation  and  ours 
agree  in  their  coefficients. 

Now  Delaunay's  expansion  has  been  recovered  by  Barton  (1967), 

Broucke  (1969)  and  Jefferys  (1970);  each  of  these  authors  report  agreement 
with  their  codes.  Therefore  we  can  presume  that  our  programs  are  opera¬ 
tional  and  our  development  is  correct. 


to  tJ  to  to  to  to 
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Table  1,  Terms  Jn 


omitted  by  Delaunay . 


Characteristics 


3 

2 

2 

1 

1 

3 

1 

1 

3 

0 

3 

3 


0 

0 

0 

0 

2 

0 

0 

2 

0 

0 

0 

0 


3 

4 

4 

5 
3 
3 

5 
3 
3 

6 
3 
3 


Arguments 


3 

2 

2 

1 

1 

1 

1 

1 

1 

0 

3 

3 


-1 

2 

_  9 

3 

3 

3 

-3 

-3 

-3 

4 
-1 
■1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

■> 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


-27 

-7 

-7 

-393 

159 

53 

-393 

159 

53 

129 

■3075 

25 


256 

32 

32 

512 

16 

256 

512 

16 

256 

130 

512 

512 


■ 

C1 

m 

m 

■ 

2 

0 

0 

6 

0 

0 

0 

ro 

35 

64 

2 

0 

0 

6 

0 

0 

2 

0 

-35 

384 

2 

0 

2 

4 

0 

0 

2 

0 

-39 

32 

2 

2 

0 

4 

0 

0 

2 

0 

-195 

128 

2 

2 

0 

4 

2 

0 

0 

0 

-15 

64 

2 

j 

0 

1 

1 

0 

0 

-261 

2  56 

> 

i 

2 

2 

1 

3 

0 

81 

16 

) 

'3 

0 

3 

1 

1 

0 

0 

27 

2  56 

2 

1 

0 

5 

1 

-1 

0 

0 
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2  56 

2 

1 

2 

3 

3 

-1 

0 

0 

81  ! 

16 

2 

3 

0 

3 

-1 

0 

0 

27 

2  56 

') 

0 

0 

6 

0 

2 

0 

0 

141 

128 

> 

2 

0 

4 

2 

1 

-2  i 

2 

0 

39 

64  ! 

2 

3 

0 

5 

5  , 

-1 

2  , 

0 

1467 

512 

2 

3 

2 

3 

1 

-1  ’ 

2 

0 

369 

32 

2  ! 

J 

0 

3 

1  1 

-1  I 

2 

0 

7011 

512  1 

9 

:  t 

5 

0 

1 

1  1 

-1 

2 

0 

2247 

512  1 

2 

1 

0  , 

5 

1  | 

-i  ! 

2 

0 

-5 

512 

2  1 

1 

2  , 

3 

1 

-i 

2 

0 

-3  1 

32  1 

2 

3 

0 

3 

1 

1 

-i 

2 

0 

-57  1 

512 

2 

i 

0  , 

5 

"I 

j 

-i 

2  1 

0 

-4401 

512 

2 

1 

2  ; 

3 

"1  , 

-i 

2 

0 

-1107  1 

32 

J  j 

0  ; 

3 

-1 

-i 

2  ’ 

0 

-4797 

512  ! 

2  i 

0  1 

5 

-1 ; 

3  , 

2  1 

0 

15 

512 

2 

1  i 

9 

1 

3 

-i 

3  I 

2  | 

0 

9 

32  ! 

2  ! 

0 

3 

-i 

» 

i  1 

9 

0 

39  1 

512 

2 

o  ! 

0  ! 

6 

0 

-2  1 

2  ! 

0 

601  1 

64 

2  i 

° 

9 

** 

4 

0  I 

-2  ! 

2 

0 

115 

4 

2  I 

2 

0  1 

4 

0 

1 

_2 

2 

0 

575 

16 

2 

0  , 

2  1 

4 

0  1 

0 

t 

0  1 

9 

45 

16  ! 

^  1 

3  i 

0  ■ 

3 

3  I 

1 

0 

0 

-2  7 

256  1 

Numerators 


Denominators 


tOWNJMNJNJNJSJIOrvJK) 
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Table  I  (continued) 


I 


* 
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3 

Tab  lt>  _J._L  To  tins  in  <t  omit  tod  by  Delaunay. 


Character  1 st  los 

Arguments 

Numerators 

Denominators 

t 

U1  |  M 

e2 

‘  1 

"2 

D1 

F1 

3 

r~ —  — 

0 

0 

4 

0 

0 

1 

0 

717 

512 

3 

(> 

(» 

4 

0 

0 

3 

0 

2115 

512 

3 

! 

0 

3 

1 

-1 

1 

0 

-33 

64 

3 

1 

0 

3 

1 

1 

1 

0 

-15 

32 

3 

» 

0 

4 

0 

1 

0 

117 

128 

3 

1 

0 

3 

1 

-1 

3 

0 

-1<>5 

8 

3 

1 

0 

3 

1 

1 

3 

0 

95 

64 

3 

1 

0 

3 

-l 

-1 

3 

0 

495 

8 

3 

1 

0 

3 

-1 

1 

3 

0 

-225 

64 

3 

0 

0 

4 

0 

_2 

3 

15325 

384 

3 

0 

0 

4 

0 

?■ 

3 

0 

5 

LO 

00 

3 

1 

0 

3 

3 

1 

1 

U 

-23 

64 

3 

1 

0 

3 

-1 

-3 

3 

0 

-7335 

64 

3 

1 

0 

3 

1 

-3 

3 

0 

2445 

64 

3 

l 

0 

3 

l 

-3 

1 

0 

-77 

32 

3 

0 

0 

4 

0 

4 

3 

0 

1029 

1024 

3 

0 

0 

4 

0 

4 

3 

0 

5 

3072 

3 

0 

4 

0 

-4 

3 

0 

177065 

3072 

3 

0 

Jj 

-4 

1 

0 

8865 

1024 

_ 

(*)  From  lioro  on  the  arguments  do  not  enter  Delaunay's  perturbation. 


-27- 


REFERENCES 

BARTON,  D.  1966,  Astronom.  J.  72_,  1281-1287. 

-  1967,  Astronom.  J.  71,  438-448. 

BROUCKE,  R.  1969,  Celestial  Mechanics  1^,  to  appear. 

BROUWER,  D.  1961,  Trans.  IAU  JJLB,  449-451. 

BROWN ,  E.  W.  1896,  An  Introductory  Treatise  on  the  Lunar  -Theory, 

Cambridge,  The  University  Press. 

DANBY,  J.  M.  A.,  DEPRIT,  A.  and  ROM,  A.  1965,  Boeing  Scientific 
Research  Labor..'  i  'es  document  Dl-82-0481. 

DELAUNAY,  C.  1860,  ThJorie  du  Mouvement  de  la  Lune,  Vol.  1, 

Memoires  de  l'Academie  des  Sciences  de  l'Institut  Imperial  de 
France,  Mallet-Bachelier,  Paris. 

DEPRIT,  A.  and  ROM,  A.  1967,  Bull.  Astron. ,  Sdrie  3,  2,  425-447. 
DEPRIT,  A.  1969,  Celestial  Mechanics  _1,  12-30. 

JEFFERYS,  W.  1970,  Celestial  Mechanics ,  to  be  published. 

MERSMAN,  W.  A.  1969,  NASA-TN-D-5459 ,  22. 

POINCARE,  H.  1909,  Lecons  de  micanique  celeste,  Vol.  2,  Part  2, 
Gauthier  Villars,  Paris. 

ROM,  A.  1970,  Celestial  Mechanics ,  to  be  published. 


